Light scattering from three-level systems: The T-matrix of a point-dipole with gain 
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We present an extension of the T-matrix approach to scattering of light by a three-level system, 
using a description based on a Master equation. More particularly, we apply our formalism to 
calculate the T-matrix of a pumped three-level atom, providing an exact and analytical expression 
describing the influence of a pump on the light scattering properties of an atomic three-level system. 
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I. INTRODUCTION 



In the description of the interaction of atomic ensem- 
bles with light, the atoms are often treated in the electric 
dipole approximation, and are therefore effectively con- 
sidered to be point-dipoles. The advantages of the 
point-dipole formalism are twofold. Firstly, the delta- 
function potential associated with point-scatterers allows 
for significant mathematical simplifications compared to 
finite-size scatterers Secondly, the point-scatterer 

formalism allows for a transparent description of many 
multiple-scattering phenomena, mimicking most of the 
associated relevant physics [2f. The light scattering 
properties of point-dipoles can be expressed by means of 
their T-matrix T which is related to the dynamic polar- 
izability a by T(cu,r',r) = —(u!/c) 2 a(uj)5(r)8(r — r'). 
If one describes the internal structure of the atoms as an 
effective two-level system Q or a damped harmonic os- 
cillator, one finds the well-known Lorentzian expression 
for the linear dynamic polarizability |4] for frequencies 
near the resonance frequency. 

In this paper, we present a extension of the point- 
dipole T-matrix to three-level atoms. Three-level atoms 
have been studied a lot in the past, especially in the 
context of lasing without inversion [5|, |6|. However, in 
this paper, we are focussing on a different application of 
three-level atoms. Our aim is to develop a transparent 
formulation of a point-dipole with gain, based on clear 
physical grounds. The introduction of a pump in the 
point-dipole model is highly attractive because it allows, 
within the framework of the T-matrix formalism, for 
a description of optically amplifying atomic systems 
@> B S TTo( | . The results presented in this paper are a 
first step towards a better qualitative and quantitative 
understanding of the presence of gain in multiple light 
scattering systems, such as atomic lattices or 
coherent backscattering experiments [l2| ■ 

Our paper is organized as follows. In Section II, the 
Master equation for an optically pumped three-level 
system is solved. The evolution of all atomic populations 
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FIG. 1: (Color online) The three-level system A. Decay from 
b to c and from c to a can occur according to the decay rates 
The and r ca , respectively. Decay from 6 to a is neglected. The 
left dashed arrows express the interaction with the pump. The 
right dashed arrows depict the interaction with the external 
probe field. 



and internal coherences can then be deduced. In particu- 
lar, we will focus on the system's steady-state regime. In 
Section III, we will connect the reduced density matrix 
of the system with its dynamic polarizability, which will 
lead us to the key result of this paper: the T-matrix of 
a point-dipole with gain. In Section IV, finally, we will 
discuss the physical meaning of our result, and elaborate 
on the dispersion and dissipation of a dipole in the 
presence of gain. 



II. MASTER EQUATION OF A PUMPED 
THREE-LEVEL SYSTEM 



We start by considering a three- level system A, shown 
in Figure [T] The experimental situation we have in 
mind is, e.g., a hydrogen-like finestructureless atom in a 
magnetic field. The c — > a transition then corresponds 
to an allowed 7r-transition, whereas level & is a higher 
lying quickly decaying level. Decay between the energy 
levels levels a, b and c can occur according to the decay 
rates r& c and r ca . The b — ► c transition either has a 
non-radiative nature, or its transition frequency is far 
detuned compared to the other relevant frequencies 
in the system. In this paper, we are interested in 
the simplest scheme possible that allows for incoherent 
pumping. Therefore, we will make the following simplify- 
ing assumptions. We assume that spontaneous emission 
from level b to level a can be neglected compared to 
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other decay processes, hence Tb a w 0. Furthermore, the 
lifetime T^ 1 of level b is chosen to be small compared 
to the lifetime T~J-. Both previous restrictions on 
the decay rates do not affect the final results of this 
paper in any qualitative way, and are only introduced 
for mathematical simplicity. Of course, quantitative 
differences will arise if extra coherences imposed by the 
pump are taken into account through a finite r^, but 
these effects are outside the scope of this paper. Two 
incident electromagnetic fields e e and e p interact with 
A. The field e p serves as a pump to establish population 
inversion in the two-level system ac. The external probe 
field e e is present in order to calculate the dynamic 
polarizability of the ac system. 

We are interested in the steady-state behavior of the 
system A. We will now derive the Master equation of 
the system, using the standard procedure presented 
in The total Hamiltonian of the system A, the 

electromagnetic fields and interactions can be written in 
the electric dipole and the rotating wave approximation 
as 

H = H a + Hp + H r + Var + V ae + Vap. (1) 

The atomic Hamiltonian Ha is given by H a b + H c , with 

H ab = fiuiS+ b S~ b , H c = hu 2 S+ c S- c , (2a) 

with S^j, S~j the ij dipole raising and lowering operators. 
The Hamiltonian of the pump is given by 

Hp = ftw p (a p al + ^, (3) 

where the operators a p and d p respectively annihilate 
and create a pump photon with frequency lo p . The 
system A is coupled via the interaction Var 

fi« --^-E(|^) 1/2 (^-4,) (4) 

k\ 

to the three-dimensional multimode electromagnetic 
field with Hamiltonian 

Hr = ^fujj k \(a k \al x + i^. (5) 

feA 

The dipole operator is denoted as d, the operators a and 
a^ are the photon annihilation and creation operators, 
and V = L 3 is the quantization volume. The external 
field e e has a frequency u>, and interacts with the system 



A by 

v AE = Inn^s+e-^ + s-^), (6) 

where the interaction strength is given by the Rabi 
frequency fl e = —d ac ■ e e /h, with d ac the ac transition 
dipole moment. 

The interaction Hamiltonian Vap, finally, denotes the 
coupling of the atom to the pump: 

Vap = 9^t b + S- b ){a p + al), (7) 

with g expressing the coupling strength between the 
pump and the system fl3j |. The pump field shifts the 
energies of the levels a and b. This shift, which can be sig- 
nificant for large pump intensities, can be taken into ac- 
count (see, e.g., [IJ]) by considering the "dressed states" 
or eigenstates of the Hamiltonian Hp, = H a b + Hp + VAP, 
given by: 

|1(A0) = sin(0)|a, N + 1) + cos(0)|b, N), (8a) 
\2(N)) = cos(0)|<z, N + l)- sin(9)\b, N), (8b) 

graphically represented in Figure [U Since the energy of 
level c is unaffected by the pump, we use the notation 
\c(N)) = |c, N). The number of pump photons is then 
given by N, and 



tan26»EE--^, < 20 < n, (9) 
d p 

where we used 

Hflp = —d a b ■ £p, S p =UJ p -LUi, (10) 

with d ao the ab transition dipole moment. We assume for 
now that the distribution of pump photons is relatively 
narrow around a large average value ((N) ^> A7V 3> 1), 
for which 2gy/ (N) — Ml p holds (see, e.g., |13|). 

The dynamics of the system A can be expressed in 
terms of the Master equation for the reduced density 
matrix a. The Master equation is in the usual Born- 
Markov approximation, written in the Lindblad form, 
given by [lE EE El 



= C nd a + C d a. (11) 
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|c, N+l> 



|b,N> 



ho>, 



|c, N> 



hco. 



la, N> 



|b, N-l> 

(a)The three-level system abc. 

|1(N)> 



a, N+l> 



tlco„ 




|2(N)> 



|1(N-1)> 



The explicit evolution equations for all reduced density 
matrix elements are found by expanding equation (fTT|) in 
the basis {|l(iV)>, \2(N')), \c, N")}, where N, N' and N" 
are not necessarily equal. However, we will show in the 
next section that only an expansion in the basis {|l(iV)), 
\2(N)), \c,N + l)} will be required to determine the sys- 
tem's T-matrix, and we will therefore restrict ourselves 
to this particular basis. If we introduce the following no- 
tations for typographical simplicity 



N,M 



E 



N 



(i(N)\a\j(M)), (14a) 

(i(N)\H D + H c \i(N)}, (14b) 
i,j e {l,2,c}, 



we can focus on the evolution equations obtained by 
an expansion of (fTTj) . Firstly, the evolution of the 
population of the dressed state |l(iV)) is given by 



■ N,N .&e . /fl w \ - 



,iwtN+l,N 



tN,N+l 



a 



l,c 



N,N 



T ca sm 2 (0)a^ N+1 - T bc cos 2 (0)<i 



sin(0) cos(8) (a{ 



N.N 



N,N 

r 2,i 



(15) 



|b, N-l> 



|2(N-1)> 



with similar expressions for the populations a 2 2* an d 



(Jc.^ ■ Secondly, coherences such as a^ 2 N evolve as 



(b)The dressed states |1(JV)> and |2(JV)). 

FIG. 2: (Color online) In Figure (a), levels \a,N), \b,N) and 
\c,N) are shown for different number of pump photons N. 
The energy differences between energy levels are also indi- 
cated. In Figure (b), the dressed states |1(7V)) and \2(N)) 
are schematically shown. The energy difference between both 
states is given by HQ = + Qp, which can easily be de- 

duced from definitions list. 



■ N,N 



■ ' e / iujt ■ in\ N+1,N —iut r n\ N,N+1 



+ T ca sin(0) cos(0)ct, 



N+l,N+l 
c.c 



be ( N,N 



1.2 



sin(0) cos(#)er. 



N,N 



sm(9) cos(9)(j^ \ N 



N 
1 1 



(16) 



The non-dissipative part of the Lindblad operator can 
be written as 

C nd a ee -~[&d + H c + V AE , a], (12) 
while the dissipative part is given by 

Ld<J = — (^S^ c S ac (T + &Sa C S a< ^J + T ca S ac crS^ c 

- V -Y + ts+X c )) + r bc s- c as+. (13) 



with analogous expressions for all other coherences, 
which are derived in the Appendix. Expanding the Mas- 
ter equation allows us to describe the dynamics of our sys- 
tem, and in particular to look at its the steady-state be- 
havior. If we examine the Master equation more closely, 
then two types of evolutions can be distinguished. First, 
the coherences are driven by the external probe field; 
therefore their evolution can be separated in a quickly os- 
cillating character (evolving typically at u>) and a slowly 
decaying envelope (evolving typically at T^ 1 ). Hence, a 
time-independent regime or "steady-state" for a(t) can 
only be obtained for the envelope of the coherences. 
The other matrix elements (the populations) however, 
do not exhibit such quick oscillatory behavior and decay 
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to their steady-state value without any persistent oscilla- 
tions. The steady-state value of the level populations 7Tj , 
i £ {a,6, c}, for example, is given in the small external 
probe field limit by 



^st _ \ ^ ^ 

77 c =2^, 



N,N 
c,c 



N 



w 



1 be st 

r b ~ 1 ' ™yi , oT c 

1 an, 



1 + W(l + 2W) 



1 - 7T 6 - 7T , 



(17a) 



where we defined the dimensionless pumping parameter 



Or 



rf, c r c 



< W < +oo, 



(18) 



of the ab transition, we can deduce that in steady-state 
(defining u>2 — w = S): 



1 - W 



(21) 

The (time-)averaged atomic ac dipole operator is 
related to the dynamic polarizability a of the ac system 
by (see, e.g., Q) 



d) — eoRe a(u>) ■ e e e 



(22) 



such that population inversion occurs in the ac system 
for W > 1. Now we can describe the dynamics of the 
system A, we will focus on the connection between the 
steady-state solution of the Master equation and the 
scattering properties of the system for radiation near 
the c — > a resonance. 



III. DERIVATION OF THE T-MATRIX FOR A 
DIPOLE WITH GAIN 

The Master equation allows us to express expectation 
values of atomic operators in terms of reduced density 
matrix elements. In the case of a single atom, the 
expectation value can be interpreted as a time-average 
due to the statistical character of the Master equation 
[HI]. The (time-)average of the operator d associated 
with the ac dipole transition, for example, is 



d) = {d-) + {d^ 

--d ac (Tr(aS~ c )+Tr(aS+ c )), (19) 



and 



Tr(as- C ) =J2H (m*s- c \m) 

N i=l,2,c 

= E(<i Ar " lsin W+< 2 Ar " cos W)' ( 2 °) 



from which we can see that indeed only the reduced 
density matrix elements given by an expansion of (jlip 
in the basis {|l(iV)}, \2(N)}, \c,N + 1}} appear in the 
expression for the average dipole moment. Furthermore, 
if we assume the pump frequency to be on resonance 



Both equivalent expressions ([19]) and ([22]) allow us to 
connect the polarizability a with the reduced density 
matrix, yielding 



a(Lu) - Bp = a nr (a„ r ■ e P x 

\ ) e ^ acy ac ^ 1 + W 

1 



-S + i^{l + W)- 



202 



(23) 



26+iT ca (l + W) 

which simplifies for small external probe fields to 



_ l-Wl 

-a 



LU 2 



1 + W 2 uj - lo 2 + + w) ' 



(24) 



The static polarizability is given by 



a = ctofj, <g) /x, 



a 



(25) 



where "(8)" denotes the tensor product of two vectors, 
and /x is the unit vector parallel to d ac . 

It is important to note that the same expression (]23|) 
is obtained if the dressing of the levels a and b is omitted 
and the optical Bloch equations are used instead of the 
Master equation. The incoherent pumping mechanism 
only appears in the polarizability as a parameter W, 
without causing any detuning effects. 

The scattering properties of a point-dipole - or, more 
generally, any scattering object - can be expressed by 
its T-matrix T (see, e.g., [l|). The T- matrix is closely 
connected to the dynamic polarizability of the scatterer. 
For a point-dipole located at r = 0, both scattering 
quantities are related by 



{r\T(uj)\r') = -(— ) a(w)5(r)5(r - r') 
— t(uj)n (§5 fi5(r)S(r — r'), 



(26) 
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with as matrix element t(ui): 



t(w) 



/w 2 \ 2 l - W 1 



uJ 2 



U-U 2 + + W)- 2Lj2 _ 2cJ +| ca(1 + W) 



(27) 



which is nonlinear in the incident probe field (through 
f2 e ). Both delta functions appearing in (|26|) express the 
local character of the scatterer; the anisotropy of the 
T-matrix is clearly due to the preferential orientation 
induced by the transition dipole moment. In the absence 
of pumping (W = 0) and for small external probe 
fields, we recover the expression for the linear dynamic 
polarizability of a two-level atom, which satisfies the 
optical theorem [2] expressing energy conservation: 



-Im 



t(co) 



uj 2 /c 



()7T 



(28) 



For the optical theorem to hold, the static polarizability 
must satisfy 



()7T 



a 



(w 2 /c) 3 LU2 



(29) 



which, if substituted in (127|) . yields the final expression 
for the T-matrix element of a point-dipole with gain for 
arbitrary pump and external probe field intensity. In the 
limit for small external probe fields, (|27|) then reduces to 



Re[t],Im[t] 
1 




(b)W = 3 

FIG. 3: (Color online) The absorption (solid line), de- 
scribed by Im[t((i))], and the dispersion (dashed line), given 
by Re[f(tj)] for the pumped three-level atom we consider, as a 
function of the detuning u — u>2 in units of T ca (a) without a 
pump (W = 0) (b) for a pumping intensity in the population 
inversion regime (W = 3). The graphs are scaled such that 
Im[t(tJ2)] = 1 in Figure (a). 



t(u>) 



3tt 



1 - W 



W 



(30) 



Expressions (f2"o| and (f3"0|) (or more generally, (f2"T)0 are 
the key results of this paper. In the next section, we will 
elaborate on their properties and physical consequences. 



IV. PROPERTIES OF THE T-MATRIX FOR A 
DIPOLE WITH GAIN 

The T-matrix (1261) fully expresses the scattering prop- 
erties of a point-dipole with gain and satisfies the 
Kramers-Kronig relations [191 ] , which can be easily veri- 
fied. 

The pump influences the T-matrix in a clear physical 
way. The most obvious effect of the pump is to induce 
changes of the (time-)averaged populations of the levels 
a and c, which leads to the multiplication of the dynamic 
polarizability with a factor (1— W)/(l+W). Remarkably, 



the multiplication factor becomes negative if population 
inversion is present in the ac system (W > I), expressing 
the fact that not only the imaginary, but also the real 
part of the dynamic polarizability is drastically changed 
by pumping. In other words, not only does absorption 
change into gain, but the dispersion relation also changes, 
as is shown in Figure [3] Besides changing the sign of the 
dynamic polarizability, the gain also effectively broadens 
the a state by a factor (1 + W), which is equivalent with 
a decrease of the Q-ia,ctor of the c — > a resonance. 

Secondly, it is intuitively clear that the optical the- 
orem does not hold any longer as a nonzero pump is 
applied. We will now show that this is indeed the case. 
If an external probe field with polarization vector e is 
incident on the point-dipole we consider here, then the 
dipole scattering cross-section a sca and extinction cross- 
section <r ext are given by 
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■40 



W 



FIG. 4: (Color online) The albedo a of the ac system as a 
function of the reduced pumping parameter W . 



&ext — 



67T 

uj/c 



l*MI 5 



Im 



t(w) 



(31a) 
(31b) 



Both scattering cross-sections depend (through t(ui)) on 
the applied pumping intensity. The albedo of the ac 
system is defined as the ratio of the elastic scattering 
cross-section and the extinction cross-section, which can 
be written for small incident probe fields as 



equal (all scattering is then elastic) in the absence of a 
pump and in the small external probe field limit. As soon 
as a pump is applied, or for larger external probe fields 
(20| . some of the light is scattered inelastically. 

Finally, we see from expressions (f3Tj) that the optical 
theorem is also satisfied for the nontrivial values W = 1 
and W +oo. The latter two pumping values cor- 
respond to pumping intensities for which the T-matrix 
vanishes (therefore the optical theorem holds), but the 
reason why is clearly different for both cases: for W = 1, 
populations are, on average, equally distributed among 
levels a and c, preventing the building up of an average 
scattered field; for W — * +oo, on the other hand, the 
linewidth broadening of the c — > a transition induced by 
the pump inhibits scattering. We stress that the absence 
of light scattering for W = 1 has to be interpreted in a 
statistical sense: the T-matrix is deduced from the (sta- 
tistical) Master equation (fTTj) , and is therefore a time- 
averaged property of the system. In other words, for 
W = 1, the pumped point-dipole does not scatter light 
on average. Furthermore, the transparency in our system 
is caused by a fully incoherent pumping scheme, contrary 
to, e.g., electromagnetically induced transparency [2lj . 

Finally, we note that the strong dependence of the 
albedo on the applied pumping intensity could lead to 
interes ting exp erimental work on, e.g., coherent backscat- 
tering [22], |23j , since the width of the coherent backscat- 
tering cone is closely related to the albedo of the scatter- 
ers in the multiple scattering medium considered. 



0~ext 



1 - W 

(i + wy 



(32) 



where (1 — a) is the fraction of the incident probe field 
which is taken away from the incident beam but not 
transformed into scattered light. If the optical theorem 
is satisfied, the albedo is one, as can be seen from ex- 
pression l|28p. Therefore, if the optical theorem does not 
hold any longer, we expect a to be smaller than unity. 
To show that this is indeed true, we plotted the albedo 
(|32p in Figure [¥] as a function of the reduced pumping 
parameter W. As soon as the dipole has internal popu- 
lation inversion (W > 1), the albedo becomes negative, 
which indicates that the point-dipole then has a negative 
extinction cross-section, as we expect. Obviously, the 
scattering cross-section is - by definition - always posi- 
tive and the effect of the pump is manifested purely as 
a decreasing of the scattered field for increasing pump. 
Furthermore, the extinction cross-section is, in absolute 
value, always larger than the scattering cross-section for 
nonzero pump. In other words, some of the incident 
probe field is taken away from the incident beam, but 
not scattered elastically. The presence of inelastic scat- 
tering is no surprise, since the total intensity emitted by 
the dipole is proportional to the population of the up- 
per state c, whereas the coherent intensity emitted by 
the dipole is proportional to the amount of coherence be- 
tween a and c (see, e.g., [l3|). Both intensities are only 



V. SUMMARY 

The aim of this paper was to find the light scatter- 
ing properties of a pumped point-dipole, modelled as a 
three-level system. The resulting T-matrix (|30[) (or more 
generally, (|27p ) is surprisingly intuitive. The influence of 
the pump can be characterized by a single dimensionlcss 
parameter W, expressing both a decrease of the Q-factor 
and a decrease of the T-matrix. Physically, the presence 
of gain not only causes the dipole to scatter partially 
inelastically, but also induces important changes in its 
dispersion and dissipation. 
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APPENDIX A: EXPANSION OF THE MASTER 
EQUATION 

The aim of this appendix is to show the expansion of 
the Master equation CCE]) in the basis {|1(JV)>, |2(JV)), 
c, N + 1)}. The reduced density matrix elements which 
are diagonal in the atomic states evolve as 



while the elements which are off-diagonal in the atomic 
states evolve as 



■N,N Me . N+1,N iut N,N+1 -iut\ 

o-i.i = -*Y sin ^H CTc ' 1 6 ~ Gi - c 6 ) 

-n ■ 2//i\ N+l.N+l T-i 2//i\ N,N 

+ r ca sm {B)(T C J ' + -r bc cos {0)u l l 

, ^ bc ■ ia\ la\( N.N , N,N\ / Ai x 

+ — sin(0) cos(0) (ct 1i2 + a h2 j, (Al) 



■ N,N Me n , a \ ( N+1,N Jut „N,N+l-iut\ 

a 2,2 = -t-^ cos{0) [a c2 e - <r 2c e J 
+ T ca cos 2 (0)<+ 1 ' Ar+1 - T bc sm 2 (0)cr?f 
+ ^sHe)cos(9)(a?f + a?f), (A2) 



. N+1,N+1 



.a 

— i — sin 
2 

.n, 

— I — cos 
2 



• f/a'i ^ n,n+i 
m(6>) ^ 

(») ( 



„ iu/t N+l,Njut 

e - cr -, e 



JV,JV+1 -iut JV+1,JV »n,t 



c,2 



JV+1,W+1 



r bf W(f?Xi 



t-i • //n m\t N+l.N+l , N+1,N+1\ 

T bc sm(d) cos(6»)(cr 1 2 



'2,1 



r 6c sm (0)cr 2 2 



(A3) 



. N.N 
°1,2 



Qe 



(sinW< 2 +1 



N e luJt - cos(9)a?f +1 e- lult ) + T ca sin(0) cos(9)a, 



N+1,N+1 
c.c 



— (<7 12 - sm(6») cos(6»)ct 2 2 - sin(0) cos^ja^i ) - ~(E? - Eg)a h2 , 



TV, TV 



N.N \ 



? N\ N,N 



(A4a) 



. TV./V+l _iu t 



a 



l.c 



^ (cos 2 (0)<r£f + V iwt - sm(0) cos(0)a ; 



/n\ N,N\ T ca N,N+1 

cos(9)a 12 J - —<J h ' c 



N,N+1 -iut 
2.c e 



j?N\N,N+l-iwt 
&1 )°l,c e 



(A4b) 



. N,N+1 m 
T 2,c e 



■ e ( m\ N+1N+1 in\ N.N ■ / n \ N,N\ lea N.N+l 

-*Y { cos ^ a ^c ' - cos(6i)a 2i2 - sm(0)a 12 j - — G 2,e 

i ^ bc ( r - ta\„„„ra\ n,n+i -iut ■ 2/ Q \ n,n+i -iut\ , * rr?N+i t?N\ 

+ -jj- [am(9) coa(9)a x ' c e - sin (9)<t 2 , c e ) + J^- c E 2 ) 



tt<N \ N,N+1 -iut 
h 2 ) (J 2,e e 



(A4c) 



and 



■ TV, AT / • JV,JV\* 

a 2,i = K^ia ) 



. N+1,N 
T e.l ' 



N,N+1,* 



. N+1,N , . N,N+U, 
°e.2 = (.^.c ) 



(A5a) 
(A5b) 
(A5c7 of 
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